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Abstract: We estimate the dominating frequencies contributing to the 
Casimir energy in a cavity of meta- materials mimicking de Sitter space, 
by solving the eigenvalue problem of Maxwell equations. It turns out the 
dominating frequencies are the inverse of the size of the cavity, and the 
degeneracy of these frequencies also explains our previous result on the 
unusually large Casimir energy. Our result suggests that carrying out the 
experiment in laboratory is possible theoretically. 
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1. Introduction 

The Casimir energy is one of the important predictions in quantum field theory and continues 
to be source of inspiration for theoretical as well as experimental work (TJ. It is the regularized 
difference between two energies: one is the zero point energy of the electromagnetic field in 
a finite cavity and the other is that in an infinite background. Generically, the expression of 
the Casimir energy depends on the details of the cavity including properties of its bulk and its 
boundary. The earliest work iff] on the Casimir energy and the following up studies J5J 01 |5J |6l 
E] |U [9] [TO] QT) all reported a result that the energy density of the Casimir energy is inversely 
proportional to the forth power of the typical size of the cavity. Applying this result to the 
universe, this kind of the Casimir energy can not be taken as a possible origin of dark energy, 
since 1 /L 4 is too small compared with the observed dark energy density lfl2l[T3"l if L is chosen 
to be the a typical size of universe. 

With doubt about the applicability of the previous results to de Sitter space, we carried 
out a calculation on the Casimir energy of electromagnetic field in the static patch of de 
Sitter space carefully lfT4l . We obtained a drastically different Casimir energy which is pro- 
portional to the size of the event horizon taking the same form as holographic dark energy 
|[T5l[T6l . Then basing on the recent theoretical and experimental development of metamaterials 
E7] [TBI [19] [20] ETJ |22l ESI |24l [25] |26] [27l , we proposed to design metamaterials mimicking 
de Sitter (along the line of studying cosmological phenomenon in laboratory, some other re- 
searches have been carried out including the negative phase velocity of electromagnetic wave 
in de Sitter, mimicking black hole and cosmic string by metamaterials [28 , 29, 30]). Different 
from the usual ones, the permittivity and permeability parameters both have divergent compo- 



nents tangent to the boundary. This unusual fact leads to a brand new Casimir energy in the 
cavity coinciding with the the result from the gravity side with the Planck scale is replaced by 
some microscopic scale in metamaterials. We encourage experimentalists to make such meta- 
materials and measure the predicted Casimir energy. This work will have significant implication 
to cosmology. 

Apparently, it seems difficult to carry out such an experiment, because the Casimir energy 
is expressed as a sum of all the frequencies, meaning that each frequency contributes a part 
to this energy, while metamaterials have frequency dispersion, in other words, the designed 
permittivity and permeability is effective only to frequencies in certain brand. However, this 
difficulty can be circumvented as we will uncover a fact that there is a typical frequency whose 
contribution to Casimir energy is dominating. Thus a cavity of metamaterials effective at this 
typical frequency is sufficient to mimic de Sitter space and induces the Casimir energy predicted 
by theoretical calculation. 

We shall show that the typical dominating frequency is co ~ l/L, where L is the size of the 
cavity, and the dominating angular quantum number is / ~ L/d, d is the short distance cut-off. 
These two numbers conspire to give the Casimir energy L/d 2 , the same order of our previous 
result Ifl4l . 

The rest of this paper is organized as follows. In sect. 2, we uncover the fact that the Casimir 
energy has a typical contributing frequency. We estimate the typical frequency in metamaterials 
mimicking de Sitter in sect. 3. We conclude in sect.4. 

2. General discussion on the typical frequency of Casimir energy 

Physically, it is conceivable that there is a typical contributing frequency to the Casimir energy. 
According to its definition, the Casimir energy measures the difference between vacuum en- 
ergy in a finite cavity and that in infinite background. The former is usually contributed by a 
discrete spectrum and the latter comes from a continuous one. Thus the Casimir measures the 
difference between discrete and continues ones. When frequency is large, the discrete spectrum 
approaches a continues one, and their contributions cancel with each other; for small frequen- 
cies, their contribution is also negligible. Since the contribution from very large and very small 
frequencies is tiny, there should be some intermediate scale at which the difference between 
discrete spectrum and continuous spectrum is maximum. Then the frequency at this scale is the 
typical frequency. 

As a heuristic example, we read off this typical frequency from the process of computing 
Casimir energy in the static Einstein's universe. In this case, the Casimir energy is given by 

ED 

E c = J-£n 3 -^ f^dco, (1) 

where the discrete spectrum consists of (0 = n/ao with degeneracy n 2 , and «o is the radius of 
Einstein universe. Reparameterizing CO by t / do, then using Abel-Plana formula 



£*■(»)- f 'dtF(t) = -F(0) + i r 'dt 
„=o Jo 2 Jo 

eq.fl} is equal to 
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(2) 
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E c =— dt-. . (3) 
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This result is very interesting because it appears that Casimir energy is contributed by frequen- 
cies satisfying a distribution similar to blackbody. We read off from this formula that typical 



frequency is at t ps 1 /2k corresponding to co ps 1/ (27Tao). However, in the general three dimen- 
sional cases of physical interest, the discussion should take into account complications due to 
at least two reasons. The first one is that it is hard to find out the exact expression of the discrete 
frequencies of electromagnetic wave in a finite cavity which are eigenvalues of a complicated 
partial differential equation deduced from Maxwell equations. The second is that in three di- 
mensions, one usually has two integers (n,l) to denote the discrete frequencies where n is radial 
quantum number and / is angular quantum number, thus Abel-Plana formula should be utilized 
repeatedly. Although it is difficult to derive an exact result, estimations are always possible 
under some reasonable assumptions. In the next section, we will try to solve the eigenvalue 
problem of Maxwell equations in the metamaterials designed to mimic de Sitter and estimate 
the typical frequency of the Casimir energy. Our estimation is based on an assumption that the 
radial quantum number n is frozen then the Casimir energy comes mostly from frequencies 
with smallest n. This is because in a finite cavity a frequency usually grows linearly with n, 
thus its contribution is suppressed exponentially by the blackbody factor. As a check of this 
assumption, the large « behavior of the eigenvalue will be given. It is indeed a linear function 
of n. 



3. typical frequency of Casimir energy in metamaterial mimicking de Sitter 

The metamaterials mimicking de Sitter space is designed with the following permittivity and 
permeability lfT4l 

£f . r= rr=L 2^mi^^ £ ee = ee = ^0_ E 99 = ^ = 1 . (4 ) 
cos(r/L) cos(r/L) cos(r/L)sin0 

where (r, 9, <p) denote the spherical coordinates. In terms of the Cartesian coordinates 

e ij = H iJ = —^JTt(S U {% sin 2 (r/L) - 1) (5) 
cos(r/L) r L r 1 

The event horizon at f = nL/2 now becomes the boundary of a cavity of metamaterials. 
The Maxwell equations in inhomogeneous medium are 

V,£' = 0, V,#<=0, (6) 

d t E l -^-djH k = Q, d t H l + ?--djE k = 0, (7) 

where is the optical metric related to permittivity and permeability through e' J = = 
ypffj, y = det(^ 7 ) and V,- denotes the covariant derivative with respect to and all the 
indices are raised and lowered by To keep the realness of the frequency of electromagnetic 
wave and the finiteness of energy, the following boundary conditions are imposed 

Ee\r=L-d — E<p\ r= L-d = 0. (8) 

These boundary conditions are acceptable physically, since the the photons emitted from the 
center of de Sitter space will travel an infinite amount of time to arrive at the horizon or they 
can never reach there as seen by any static observer. 

To solve Maxwell equations we adopt Newman-Penrose formalism [32]. That is to use four 
null vectors reexpress the Maxwell tensor F^ v as 



Ffiv = 2[<j>i{n^l v] +m^m* v] ) + ^l^m v] + <j>om*^n v] }+ c.c, 



(9) 



where "[ ]" denotes the antisymmetrization, and "c.c" means the complex conjugate. The con- 
vention about s and the null vectors is given by 

00 = F^m v , 0! = if)xv(' M « V +m*V) , 2 = F^ v m*»n v , (10) 



with 



= (- 1 TT ^,l,0,0), m» = -L(0,0,1,- ' 



l-r 2 /L 2 ' ' ' ' ' V2r ' ' ' sin0 ; 

» M = (J, -^-^,0,0), m** = -L(0,0,1,— U (11) 

To solve Maxwell equations conveniently, we have adopted a coordinate system different from 
that appearing in dl}, but we will transform back to the old coordinates when finding the typical 
frequency. Then after some standard steps 11331 . we obtain 



</. 1=e --'y/"(0,(p^(r), (12) 
where Yl"(6,(p) is the spherical harmonic function satisfies the following equation 

dl 



sin 6 sin 



(13) 



and R(r) satisfies 



2 2 2 2 2 

[r 2 (l - ^)d 2 + 4r(l - ^ H + T— ^7 - + 2- /(/ + l)]*(r) = 0. (14) 

This equation possesses two independent solutions and the one of physical interest is 

R(r) = r 2 {\ + -)-<- 1 (1 + -)-' wL l 2 (\ - -) iwL l 2 F{l + 1,1 + l + iwL,2l + 2, —), (15) 

it is regular at r = and goes back to fiat space result ji{wr)/r when L — >• oo, The solution of 
00 and 02 can a ls° be found, but to find out a> solving 0j is enough. 

Recall that in an empty spherical cavity the two independent electromagnetic modes are 
TE and TM corresponding to E r = and H r =0 respectively. Here the situation is similar. In 
the TE modes, E r = the electrowave is transverse. From the definition of <j>s eq.dTOl). we 
read 0i = — \ (E r + iH r ) = — l -H r . Combining the following Maxwell equations and boundary 
conditions (© 

iwH r = 1 {d e E<p - d<pE ), E e \ r= L-d=E (p \ r=L _ d = Q, (16) 
r z sin0 r r r 

We deduce that 

h\r=L-d=0. (17) 

Since d <C L, L — d ~ L, the boundary conditions can be imposed on the r — > L behavior of 0i . 
When r — > L, the radial part of 0] has the following asymptotic form 



Then eq.dTTb requires that 



F(-iwL) \-r/L iwL/2 
S{ T{l + \-i W L) { \+ r /L } nr=L - d U ' UyJ 

In the coordinate system eq.(5), the cut-off d is imposed on the physical radial coordinate f ©, 
the cut-off on r is therefore of 2 / (2L), the above condition amounts to 



T(-iwL) A-r/L iwL/2 
M T(l+l-iwL) { TT7jL> >\r=L- d y 2 L-0 (20) 



which determines value of CO. For small /, we can solve eq.d20l) to pick out the lowest CO when 
ln(4L 2 /d 2 ) > 1 (This is guaranteed by the fact that in usual metamaterials d is nanometer, and 
L is 1cm) and the other corresponds to very large n. The results are exhibited below. 

1) / = 0. We find that (EUJ leads to 

sin[f H AL 2 /d 2 )] = => co = Lln ^ 2/d2y n = 1,2- • • . (21) 

2) / = 1. Then d20j implies that 

COL = tan( — ln(4L 2 /W 2 )). (22) 

Since \n(AL 2 /d 2 ) > 1, the lowest CO « 27r/Lln(4L 2 /c/ 2 ). 

Based on the above results, we infer that the Casimir energy from the frequencies correspond- 
ing to small I is of the order 1 / L\n(4L 2 / d 2 ) which is much smaller than L/d 2 predicted in fl4l . 
So the typical frequency cannot be around 1 /Lln(4L 2 /af 2 ). To estimate the typical frequency, 
we observe that for / 1 the Stirling formula can be used to reexpress eq.d20b as 

Re[r(-/a)L)(£///2L) ,fl)i ] = 0. (23) 

We see that a critical I denoted by l c emerges at l c = 2L/d (Indeed l c ^ 1, so our estimation 
is reasonable). For I ^> l c , the term (dl/2L)"° L becomes highly oscillating, and the CO satisfy- 
ing d23l approaches a continuous distribution whose effect is canceled by that from the infinite 
background. For I <C / c , this case is just what we discussed before, their contribution is sublead- 
ing. Therefore the dominating contribution to the Casimir energy can only come from I ~ / c . 
When / l c , the corresponding frequency is around 1/L since in this case the constant appear- 
ing d23l is of order 1 . Go one step further, we estimate the contribution to the Casimir energy 
from frequencies corresponding to I ~ l c . For l c /2 < I < 3/ c /2, for each / the degeneracy is 
2/ + 1, 

3/ c /2 

E c ~ £ (21 + 1)1/1 ~ l 2 /L = L/d 2 . (24) 

k/2 

This result is the same order as we obtained in lfT4l by a different method. It is a strong support 
to our estimation about the typical frequency. 

In the TM modes, the magnetic wave is transverse 

H,=0, fr = ~E r . (25) 
Combined with the Gauss law and the boundary conditions (O 

[d r (r 2 sin 0E r ) + ^^d e E e + ±- d^} \ r=L _ d 2 /2L = (26) 

1-77 (I-77) sin0 



We deduce that 



dr{r^)\ r=L _ d 2 l2L = d r (r 2 E r )\ r=L _ d2/2L = 0. (27) 
After some simplification, this condition is transformed to 
W r(-iwL) , l-r/L iwL/2 

lm{ T(l + l-i W L) i TT7jL ) ' )\r=L-<P/2L-0. (28) 

Then the process of finding frequency CO satisfying above condition is the same as before and 
we list the results as follows 

1) / = the frequency is given by 

,coh , n. , ... (2n + 

cos[— ln(4L 2 /d 2 )} = 0^co = ^^j— ,„ = 0,1-. (29) 



2)1 = 1. 



coL=-cot(^\n{4L 2 /d 2 )). (30) 



For ln(2L/d) ^> 1, the lowest frequency is 

co^7t/Lln(4L 2 /d 2 ). (31) 

3) The typical frequency can be read from 

ha\T{-i(0L){dl/2L) imL ] = 0. (32) 

Since the structure of eq.d32l) and eq.(l23l is similar, the typical frequency is still of the order 
1/L with the critical l c given by 2L/d. 

We emphasize that to estimate the typical frequency we have assumed that the contribution 
from co with large radial quantum number n is suppressed exponentially. As a check of this 
assumption, we present the expression of CO in large n limit. Both TE and TM modes have the 
following frequencies 

m "LW™ 2 /diy (33) 

Thus for n is large, co grow with n linearly and their effects will be suppressed by the black 
body factor. 

To end this section, we propose that the measurable quantity for such experiments would be 
the Casimir force. From eq. (l24b . we get 

Fc--^_- (34) 

It is remarkable that this attractive force is relatively large compare with the usual Casimir force. 
Thus, it is more easily to measure the new Casimir force for a cavity with small cutoff. We also 
notice that the unusually large Casimir force is related to the permittivity and the permeability 
atf = f {L-d) © 

£ = M~§- (35) 
d 

where d is some microcosmic cutoff to keep the permittivity and the permeability finite but 
relatively large. 



4. Conclusion 

We estimate the typical frequency of the vacuum fluctuations in metamaterials mimicking de 
Sitter and find it proportional to l/L, the size of the cavity. Assuming d be 1 nanometer, and L 
be 1cm, then the typical wavelength is about 1cm. 

With our estimation of the typical frequency and the typical angular quantum number, we 
also have an intuitive understanding of our Casimir energy formula. 

We hope that one day the metamaterials suggested by us can be made with appropriate size 
and effective for the corresponding typical frequency, then the predicted brand new Casimir 
force can be measured, this experiment is important for study cosmology in laboratory. 
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